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Atoms in highly excited Rydberg states exhibit remarkable properties and constitute a powerful tool for studying quantum 
phenomena in strongly interacting many -particle systems. We investigate alkali atoms that are held in a ring lattice and excited to 
Rydberg states. The system constitutes an ideal model system to study thermalization of a coherently evolving quantum many- 
particle system in the absence of a thermal bath. Moreover, it offers exciting perspective to create entangled many-body quantum 
states which can serve as a resource for the generation of single photons. 



1 Introduction 

During the last decades, ultra cold atomic physics has experi- 
enced a tremendous boost initiated by the invention of laser 
cooling. In this course, and in particular with the exper- 
imental achievement of a Bose-Einstein Condensate, it was 
quickly realized that ultra cold atomic systems have the po- 
tential to explore phenomena across many areas of modern 
physics. This is rooted in the fact that these systems are ex- 
tremely versatile for their interaction properties as well as their 
confinement can be tailored by optical and magnetic fields to 
great precision^. As a consequence, many model systems that 
served our understanding of many-body phenomena in con- 
densed matter physics became suddenly realizable and could 
be studied experimentally. Prime examples, to name only a 
few, are the experimental investigation of the Mott-Insulator 
transition^, the Beresinskii-Konsterlitz-Thouless phase transi- 
tion in a two-dimensional bose gas^, or the recently observed 
pinning quantum phase transition in a Luttinger liquid^ . These 
and further experiments have had a huge impact on the ad- 
vancement of modern physics across area boundaries and have 
stimulated further work in condensed matter physics, quantum 
optics, quantum information, ultra cold chemistry and statisti- 
cal mechanics. 

1.1 Rydberg atoms - Excited states with exciting perspec- 
tives 



At present the majority of ultracold atoms experiments is car- 
ried out with ground state atoms. Very recently, however, there 
is a growing initiative towards exploiting the unique properties 
of atoms in highly excited states. These so-called Rydberg 
atoms are blessed with remarkable properties'^. Firstly, albeit 
highly excited, their lifetime can reach tens or even hundreds 
of microseconds. This is to be contrasted with the typical life- 
time of the first excited state of alkali metal atoms which is on 
the order of several tens of nanoseconds. Secondly, atoms in 
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Rydberg states interact strongly over large distances - the in- 
teraction strength can reach several tens of MHz at a distance 
of the order of a micrometer ^ . Clearly, the involved timescales 
and interaction strengths are vastly different from those en- 
countered in 'traditional' ultra cold atoms setups where the 
interatomic interaction is well- approximated by a contact po- 
tential and interaction energies are of a few kHz. These hugely 
different properties illustrate that many -body physics with Ry- 
dberg atoms takes place in a completely different parameter 
regime, which is usually referred to as 'frozen gas'^^. Here 
the quantum dynamics occurs in the internal atomic degrees of 
freedom - the electrons - while the atoms itself can be regarded 
as being fixed in space. 

Experimental studies of this excitation dynamics have been 
performed by a number of experimental groups; predomi- 
nantly in setups where Rydberg states were excited from an ul- 
tra cold gas of alkali metal atoms. All these experiments show 
a dramatic reduction of the fraction of excited atoms once the 
atomic density or the interaction strength surpassed a certain 
value This is a manifestation of a central hallmark in 
many-particle Rydberg physics - the 'Rydberg blockade '1^^^. 
This blockade originates from the strong interaction among 
Rydberg atoms which inhibits the simultaneous laser excita- 
tion of two nearby atoms to Rydberg states due to the induced 
level shifts. It is this effect which is responsible for the rich dy- 
namical behavior and collective character of ensembles 
of Rydberg excited atoms. In very recent experiments the Ry- 
dberg blockade has been experimentally demonstrated for two 
atoms which were held in separate optical trapsl^^^. More- 
over, the implementation of coherent quantum operations be- 
tween two such atoms has been successfully shown This 
supported, in an impressive fashion, the feasibility of vari- 
ous theoretical proposals that use Rydberg atoms for the im- 
plementation of quantum information processing schemes^, 
the study of the evolution of strongly correlated quantum sys- 
temsl ^^tt25 | coherent quantum state preparation^^^^ and quan- 
tum simulation 

In a very recent development, Rydberg states are about to 
find their way also into conventional ultra cold atoms experi- 
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ments as their properties can be used to create long-ranged in- 
teractions^^. In order to overcome the above-mentioned mis- 
match between lifetime and interaction energy Rydberg states 
are merely weakly admixed by a far off-resonant laser The 
thereby created interatomic interactions will potentially per- 
mit the observation of exotic quantum phases such as super- 
sohdsEl. 

1.2 Outline 

In this work we will review recent and new theoretical work 
carried out on so-called 'Rydberg rings'. In the system that 
is central to this work, ultra cold atoms are confined to tight 
traps which are arranged on a ring as shown in Fig. [T] Upon 




Fig. 1 In our setup the atoms are held in L tightly confining traps 
which are regularly placed on a ring with radius R. The typical 
distance between the traps is a few micrometers. 

laser excitation the atoms interact strongly and the subse- 
quent real time dynamics as well as the static properties will 
be thoroughly discussed in this paper. The system consti- 
tutes a paradigm example for the versatility of the Rydberg 
atoms. We will demonstrate that the system is ideally suited 
to study the transition of a clos ed qu antum many-particle sys- 
tem into 'thermal equilibrium'ISISl, Understanding the gen- 
eral mechanism that underlies the the thermalization of closed 
many-particle systems is currently a very active direction of 
researching. Moreover, we will show t hat it allows to cre- 
ate entangled many particle states^^'^^'^^' which can serve as 
resource for the creation of single-photon light sources——. 

2 The Rydberg Ring 
2.1 The ring lattice 

Our system is formed by bosonic atoms confined to a ring 
shaped one-dimensional lattice. The L lattice sites are well- 



separated - with inter- site spacing a - so that the low-lying 
quantum states of each site can be approximated by eigen- 
states of a local harmonic oscillator potential. The lattice spac- 
ing is of the order of several micrometers which in principle 
allows for single site resolution with moderate experimental 
effort. Such lattice can be realized by using deep large spac- 
ing optical or magnetic lattices Throughout this work 
we will focus on a regime in which on each site the atoms pop- 
ulate the harmonic oscillator ground state, whose spatial width 
a is much smaller than the lattice spacing a. In fact, we will 
assume that the atoms are infinitely localized, i.e., a/a 
Moreover, we will consider that each site contains an identical 
number of atoms which we denote by A/q. This very idealized 
situation can only be approximately achieved in experiment. 
A fluctuating particle number and a finite width of the local- 
ized atomic wave packets caused by the quantum uncertainty, 
interactions and a finite temperature can in principle be ac- 
counted for. All these effects introduce disorder in the system, 
which slightly alters its properties'^. This is a very exciting 
topic by itself and will be expanded on elsewhere. 

2.2 Level structure and interaction 

Let us now focus on the electronic level structure of the 
trapped atoms. In this work we consider Rydberg states of 
alkali metal atoms whose level structure is particularly simple 
and strongly reminiscent of the hydrogen atom^. States with 
low angular momentum, however, are - unlike in the hydrogen 
atom - energetically well separated from a highly degenerate 
manifold of levels that is formed by states with high angular 
momentum. 

In experiment Rydberg states are usually excited by a two- 
photon transition (see e.g.EEH)^ Here the ground state of 
the atom \g) is coupled via an intermediate /?-state to a Ry- 
dberg state |r). Under certain conditions which are usually 
satisfied experimentally the intermediate state can be elim- 
inated and the two-photon transition is effectively replaced 
by a single fictitious lasei^. Due to the selection rules for 
dipole radiation only s ov d Rydberg states can be excited 
this way. In this work we choose |r) to be an ^-state, for 
these states possess two peculiar properties. First, they are 
energetically well isolated. Second, atoms excited to these 
states are, to a very good approximation, interacting via the 
isotropic van-der-Waals interaction. This interaction is of the 
form yvdw(r) = C6 x |r|~^ where r is the interatomic separa- 
tion and C6 parameterizes the interaction strength. Of impor- 
tance here is the scaling of with the degree of excitation 
which goes with the eleventh power of the principal quan- 
tum number^^^. As a consequence the van-der-Waals inter- 
action, which is ubiquitous among atoms, is strongly exagger- 
ated when Rydberg states are involved. 

This interaction produces a severe alteration of the excita- 
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' van-der-Waals interaction 




Fig. 2 a: Level scheme for two interacting Rydberg atoms where 
the level structure of each atom is approximated by the two states 
\g) and |r). If the distance between the atoms is smaller than the 
blockade radius rb the doubly excited state is no longer accessible 
(see text), b: If an ensemble of atoms is confined to a spatial region 
in which the distance between all atoms is smaller than r^, the laser 
can only excite a single atom at a time. This creates a superposition 
of all possible possible states that carry a single excitation, often 
referred to as super atom. 



tion properties of atoms confined to a single site o f the lat- 
tice which is known as the Rydberg blockade^SE]^ Xo un- 
derstand this let us consider the situation depicted in Fig. 
[2^, i.e., two atoms separated by a distance r whose states 
\g) and \r) are coupled resonantly via a laser with Rabi fre- 
quency The energy of the doubly excited states |r)2 
is twice the atomic transition energy plus the interaction en- 
ergy y(r) = C6 X |r|~^. This additional energy has to be over- 
come for the simultaneous excitation of two Rydberg atoms. 
This can only happen if the line corresponding to this transi- 
tion (sketched by the blurred level in the figure) is sufficiently 
broad. For a strong laser the line width is determined by the 
Rabi frequency. This means in turn that atoms cannot be si- 
multaneously excited if their separation is smaller than 
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i.e., if the interaction induced energy shift is larger than the 
Rabi frequency. The quantity rb is called the blockade radius. 
Throughout we will assume that this radius is much larger than 
the diameter of a single lattice site. Thus, out of of A/q atoms 
located on a given site k only a single one can be excited, 
which is sketched in Fig. ^p. This means that on each site 
only two (collective) states are accessible, i.e.. 
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where 5 is the symmetrization operator. The state |/?)^ is 
sometimes called superatom^^ and consists of a symmetric 
superposition of all possible single Rydberg excitations. 

Let us now look in more detail on the excitation behavior of 
such a single blockaded site. The action of a laser with single 
atom Rabi frequency Q.^ and detuning A is governed by the 
Hamiltonian 



Ho = E 



k=\ 



^o\r)k{g\^^\r)k{r\- 



-h.c. 



Having only two possible states on each site reduces the prob- 
lem to that of a spin system where |G)^ = ||)^ and \R)j^ = W)^. 
Thus, in the restricted subspace spanned by the collective 
states ([3j, the Hamiltonian becomes 



where 



L 

^laser ^ 
k=\ 
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(k) (k) 

Cy and cl being the Pauli spin matrices. The laser coupling 
is strongly enhanced due to the collective nature of the state 
\^)k ~ It)^- The corresponding collective Rabi frequency is 
given through = Qq^/Nq. 

We are now in position to formulate the spin Hamiltonian 
that governs the excitation dynamics of the Rydberg ring. Un- 
der the assumption that the atoms are strongly localized at the 
lattice sites the interaction between nearest neighbors is given 
by V = Ce/a^. Due to the strong decay of the interaction with 
distance, it is sufficient to consider only this nearest neighbor 
interaction. This is a valid approximation provided that the in- 
teraction is not so strong that the blockade radius encompasses 
next-nearest neighbor sites, i.e., > V /64. The Hamiltonian 
of the system then becomes 



with the interaction Hamiltonian 



(5) 



k=l 



This constitutes the basis for the remainder of this paper. 

Throughout we will not consider radiative decay of the Ry- 
dberg atoms. This is justified if there is a clear separation 
between the timescale at which the laser-driven dynamics is 
taking place and the atomic decay rate. In a typical exper- 
imental situation such separation is possible: The dynamics 
takes place on a timescale ranging from the 100 nanoseconds 
to a few microseconds^^ while the lifetime is on the order of 
50 to 100 microseconds. 
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3 Real time dynamics and thermalization 

In this section we will study the real time dynamics of the 
model ^ in the strongly interacting regime, i.e., when V :^Q,. 
We start from an initial state in which all atoms are in the 
ground state, that is, the state of the system is given by the 
product |init) =Y[k\^)k- ^i^^ show that the long time be- 
havior of the mean number of excitations can be understood 
through equilibrium thermodynamics of an equivalent classi- 
cal dimer model. 

3.1 Long time behavior of observables and perfect block- 
ade model 

In Fig. |3]we show the temporal evolution of the mean number 
of excited particles, N = Y,k^k and the density-density corre- 
lations, gnm = i^m^n) / {N)^ bctwccn neighboring sites and 
next-nearest neighbors. The data shown is in fact calculated 
by considering fully the 1/r^-tail of the van-der-Waals poten- 
tial and an interaction strength V = 5Q,. In this regime the 
next-nearest neighbor interaction employed in eq. ^ pro- 
vides an excellent approximation. All three displayed quanti- 
ties show similar behavior: For short times we observe a num- 
ber of large contrast oscillations whose amplitude diminish for 
longer times, so that the quantities approach a steady state, 
i.e., all of them assume a quasi time-independent value with 
only very small amplitude fluctuations. Our aim is to under- 
stand the physical origin of this steady state which, however, 
occurs only for sufficiently strong interaction^^ To this end 
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Fig. 3 Temporal evolution of the mean excitation number (N) {t) 
and the correlation functions guit) and gi3{t) (L = 20,V = 5Q). 
For short times the data shows pronounced oscillations of large 
amplitude which are followed by a steady state regime in which 
small amplitude fluctuations about a stationary value take place. 



are accessible by a time evolution. This in turn means that 
two excited atoms cannot simultaneously occupy neighbor- 
ing sites, i.e., within this model the correlation function gi2 is 
strictly zero. For resonant laser excitation (A = 0), which we 
consider throughout this section, this model has no adjustable 
parameters and by scaling time SiSt ^ Q,t the energy scale be- 
comes unity. In this regard it constitutes a fundamental model 
for the understanding of the dynamics of a strongly interacting 
one-dimensional Rydberg gas. 

3.2 Evolution in excitation number space 

The strongly interacting case is characterized by a pronounced 
collective behavior of the system. That means the 'good' de- 
grees of freedom or eigenexcitations of the Hamiltonian are 
very different from single atom excitations. Hence an under- 
standing of the long-time behavior is most likely to be ob- 
tained by looking on the system from an angle that puts less 
emphasis on single atom properties. 
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Fig. 4 a: Sketch of the system's structure in excitation number 
space. The laser connects configurations whose excitation number 
differs by one. The precise strength of the connection is contained in 
the entries of the matrices b: Transitions can only take place 

between compatible configurations, i.e., the laser can only connect 
states that can be converted into each other by the addition or 
removal of exactly one excitation. 



we employ a model which is even simpler than eq. ([5]), called 
the perfect blockade model Its essence is that the nearest- 
neighbor interaction is assumed to be infinite (V ^ oo). As a 
consequence, only the quantum states |<I>) which satisfy the 
constraint 

//int|<^>)=0, 



Our approach to the problem of understanding of the steady 
state value of the mean excitation number is to study the evolu- 
tion of the system in excitation number space. The underlying 
insight is that the configurations of the system can be charac- 
terized by two properties: Their number of excitations and the 
actual way these excitations are distributed on the ring. The 
laser can only couple configurations whose excitation number 
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differ by one and which are compatible in the sense that the 
removal/addition of exactly a single excitation converts one 
configuration into the other. This can be cast into a diagram 
which is shown in Fig. |4^. An example for a two compatible 
and two incompatible states is given in Fig. [4|). Our aim is 
to understand the dynamics along the horizontal of this dia- 
gram as this will show us directly how the excitation number 
is distributed. 

In the above representation the Hamiltonian acquires a 
block structure where each block corresponds to a different ex- 
citation number n. Since the laser is on resonance, all entries 
of each of these square blocks of dimension dim„ x dim„ are 
zero. Here, dim„ is the number of configurations that contain 
exactly n excitations. The blocks are coupled through rect- 
angular submatrices Cn,n+\ with dimension dim„ x dim„+i. 
Thus, defining the projector on the subspace of n excitations 
as 1^) (^1 the Hamiltonian can be written as 

^max 

H = \m){m\H\n){n\ 

n,m=0 



£ Cn^n^i\n){n^l\^Cl,^,\n^l){n\ 

n=0 



(6) 



where {n\H\n-\-l) = • Here Umax is the maximum num- 
ber of excitations which is given by L/2 and (L — 1 ) /2 for even 
and odd lattice size L, respectively. 

3.3 The rate equation 

We will now show that it is possible to derive an equation for 
the system's evolution in the excitation number space that pos- 
sesses a steady state solution which explains the saturation of 
the number of Rydberg atoms observed in Fig. [3] 

The mean number of Rydberg atoms that are excited after a 
time t is given by 



{N){t)='£npn{t), 

n=0 



The first term vanishes if p(0) is diagonal in the excitation 
number space, i.e., if the initial state contains a defined num- 
ber of Rydberg atoms, which we assume here. The second 
term can be rearranged such that 

dtPn = l^dxTr [{2H\n){n\H-{\n){n\ ,H^}) p(x)] . (9) 

Let us study the first term of this equation in order to see how 
it can be further simplified: Upon insertion of two complete 
basis sets and abbreviating pmk{^) = p{^) 1^)1*]^^ 



Y^TrWk) {k\ H \n) {n\H \m) (m| p(x)] = £Tr C,,„4,„Pm^ W- 
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The crucial observation is now that Cj^^C^^ ~ Kmn^m^km 
where K^n is a constant and 1^ is the dim^ -dimensional rep- 
resentation of the identity matrix that at the same time is the 
projector onto the subspace containing m excitations, so that 
Tr 1 

mpmm] — Pm- The reason for this is that the matrices 
are sparsely occupied and that the entries, due to the strong 
interaction, are completely uncorrelated. Thus, when calcu- 
lating the products Q^Cji^^, the only elements that acquire an 

appreciable size stem from a row of Q ^ being multiplied by 
itself, i.e., the corresponding column of C^^^. This happens 
only if m = ^ and only for the diagonal elements. This is il- 
lustrated in Fig. 5 where we look at the product Q ^C^^ for 
L = 20. We thus Tind 

^56 




\ 






\ 



(7) 



Fig. 5 The product C^^C^^ (L = 20) is approximately proportional 
to the identity matrix with dimension dims . 



where pn{t) is the probability to find exactly n excitations. 
This probability is given by Pn{t) = Tr[|^)(^| p{t)] with p(^) 
being the density matrix of the system. The evolution of 
the density matrix is governed by the von-Neumann equa- 
tion dtp{t) = —i[H,p{t)] which is equivalent to the integral 
equation p{t) = p(0) — //q Jx[if ,p(x)]. Inserting this inte- 
grated equation into the von-Neumann equation and abbrevi- 
ating pn{t) = Pn one finds 

dtPn = -iTr[[H,p{0)] \n){n\]- T JxTr [[//, [//,p(x)]] \n){n\]. 

Jo 

(8) 



Y,^^ Q,n^knPmk{^) -k^-i,^/?^-i(t:) + k„+i,^/?^+i(x). (10) 

km 

There are three more steps which we have to perform in or- 
der to arrive at a closed rate equation for the probabilities 
Pn{t): Firstly, we treat the two remaining terms of eq. ^ 
in the same way we have treated the first term. Secondly, we 
make use of the fact that Tr [C^ ^Cji = Tr [cj^ which 
fixes the constants K„,^+i = w„dim„+i and K„+i,^ = W/^dim^ 



* Note that Pmki^) is not a matrix element but in general (due to the large de- 
generacy of the excitation number subspaces) a dim^ x dim^-matrix. 
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where the are unknown coefficients which we assume are 
non-zero. Thirdly, we ignore the behavior of eq. ^ on very 
short time intervals, i.e., the evolution of the system on a time- 
scale shorter than which is the typical timescale given by 
the laser. Assuming that during such a coarse grained time- 
interval of length X the density matrix does not change we ar- 
rive at the following equation for the probabilities 



Pn{^)-Pn{0) 



X lUfi— 

idimnPn-l UndimnPn^l] 
-X[Undimn^l^Un-ldimn-i]pn. (H) 



A more quantitative and thorough discussion of the derivation 
of this equation is given in Ref.^"^. 

The steady state of this equation satisfies 



Pn 



steady 



and is given by 



steady _ '^^n 

dim 



(12) 



with the dimension of the total Hilbert space dim = 
Ln=o pT^^^ represents the probability distribution of 

the number of excited atoms or the full statistics of the Ryd- 
berg number count in the steady state. Remarkably it can be 
obtained without a precise knowledge of the numbers 



For the Rydberg ring the steady state solution ([12]) can be 
given explicitly: 



with 



steady . 
Pn 



dim : 



1 



dim L — n 



L — n 
n 



L — n 
,L-n \ n 



(13) 



(14) 



This allows us to extract the mean value of the number of Ry- 
dberg atoms in the steady state, which for L ^ oo evaluates 
to 



1 



1 
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^ 0.276L. 



(15) 



3.4 Numerical results 



Let us now compare our findings to the numerical solution of 
the system when the initial state is the product state |init) = 
Y\k\^)k' i-^-' /^o(0) = 1- In Fig- [6^ we show the temporal 
evolution of the probability density pn in excitation number 
space, i.e., along the horizontal axis of the diagram shown in 
Fig. |4^. For short times we observe oscillations which quickly 
diminish after ^ 5. Beyond that time the distribution be- 
comes strongly localized and maintaining its shape with slight 
oscillations on top. Here the system is localized in excitation 
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Fig. 6 Temporal evolution of the probability density pn in a system 
consisting of L = 25 sites in excitation number space. Initially all 
atoms are in the ground state, i.e. = 1- a: In the interacting case 
the system reaches eventually a state in which the probability 
density localizes in excitation number space, b: This is not the case 
in the absence of interactions. Here, the wave packet performs 
coherent oscillations with maximal amplitude. 



number space populating mainly configurations with excita- 
tion numbers close to (A^) . This relaxation occurs due to the 
quasi-random or uncorrected couplings between the different 
excitation number subspaces which constituted the basis for 
the derivation of eq. ( pT] ). This eventually washes out all phase 
coherence between the subspaces containing a different num- 
ber of excitations. In the non-interacting case this phase co- 
herence is maintained throughout so that the population shows 
(Rabi-)oscillations of full contrast, as shown in Fig. [6b. 




Fig. 7 Probability density pn in excitation number space for L = 25. 
The green (thin) curves are snapshots taken during the interval 
100 < ^ < 104. For these times the calculated Rydberg number 
shows the steady state shown in Figs. [3]and[6^. The dashed curve is 
obtained by taking the average over the set of snapshots. The red 
thick curve shows pn^^^^ as given by eq. |l3| . 

To compare the numerically obtained distribution to the an- 



alytically predicted one ([13]) we remove the temporal fluctua- 
tions (which are a finite- size effect) by averaging the distribu- 
tion over a time interval of a certain length. We here choose 
the interval 100 < < 104 during which the steady state is 
well-established. The result is presented in Fig. |7] The agree- 
ment is good, both distributions are peaked at the same value, 
i.e., n = 1 for L = 25. However, the numerical result has a 
systematically enhanced probability at low n. This deviation 
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is also reflected in the mean number of Rydberg atoms which 
yields (^^) numerical ^ 0.26L which is slightly smaller than the 
result ( p3] ). The reason lies in the particular choice of the ini- 
tial state. When choosing a random initial state (which might 
be difficult experimentally) one obtains in general a perfect 
agreement between the numerical data and the analytical re- 
sult as is shown in Ref.'^. 

3.5 Connection to the classical hard core dimer model 

Having found strong evidence for the fact that our coherently 
driven and closed quantum system thermalizes, it is tempting 
to ask the question whether there exists an analogous classi- 
cal model whose thermal equilibrium properties coincide with 
the properties of the quantum model for long times. In our 



mined solely by the entropy. This means that the probability 
Pn is simply proportional to the number of available configu- 
rations or microstates that contain n such excitations. 

We have to emphasize at this point that the quantum state 
of the system at any time possesses zero entropy, even in the 
steady state. This is a consequence of the fact that a pure quan- 
tum state always remains a pure state under a coherent evolu- 
tion. The entropy we referred to previously is the one related 
to the classical subsystem. This classical subsystem is charac- 
terized by a density matrix in which all the information which 
is not necessary for evaluating properties of classical observ- 
ables - such as the mean particle number or density-density 
correlations - is traced out. In the steady state this density 
matrix has maximum entropy, i.e., it is completely mixed. 



Fig. 8 Classical lattice gas of hard-core dimers whose thermal 
equilibrium properties coincide with the steady state that is assumed 
by the Rydberg ring for long times. The laser adds or removes 
dimers which cannot simultaneously occupy the same lattice site. 

system the laser adds/removes Rydberg atoms to/from the lat- 
tice. The size of the Rydberg atoms is effectively two lattice 
sites and they are also hard objects as there is strict nearest 
neighbor exclusion. This is reminiscent of a classical hard 
core dimer model as shown in Fig. [8] Let us now calculate the 
mean number of dimers occupying the lattice in this classical 
model. To this end we first assume that the dimers are not hard 
objects but that the simultaneous occupation of a site by two of 
them costs an energy V. The grand canonical partition func- 
tion of this system is then given by S(|3,/i,L) = (^+)^ + (^-)^ 



with ^4 



^± = 2 (l j ±^^J4e\^^^{l-e\^^^-^)y where 

P is the inverse temperature and jn the chemical potential. The 
mean number of dimers in equilibrium is then given by 
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4+(e-Pi'-l) 



pv- 



1 



V5, 



(16) 



which, in the limit PV oo, i.e., infinitely strong interaction, 
becomes the result that we have previously obtained for the 



number of Rydberg atoms [eq. (15)]. This is not surprising 
because eq. ([14]), i.e., the dimension of the Hilbert space of the 
perfect blockade model, is by construction the microcanonical 
partition function of the hard core dimer model or onedimen- 
sional monomer-dimer problem 

This shows that in equilibrium the probability for the sys- 
tem to reside in a subspace with n Rydberg atoms is deter- 



3.6 Beyond the perfect blockade model 

We have seen that the microcanonical ensemble very accu- 
rately predicts the long time properties of the distribution of 
Rydberg excitations. All configurations occurred with equal 
weight which is a result of the fact that we chose to con- 
sider only configurations with zero interaction energy. When 
accounting for a finite but large nearest-neighbor interaction 
strength V, a significantly larger number of configurations be- 
comes accessible during the dynamics. The population of con- 
figurations with high interaction energy will, however, be sup- 
pressed during the time-evolution since those are far detuned 
from the initial state with zero interaction. Interestingly, it 
turns out that in the long-time limit the properties of the sys- 
tem can again be understood by using statistical mechanics 
if the interaction strength surpasses a certain value which is 
Vtrans ^ 2Q.^^ . Here, for example, the mean particle number 



is given by eq. ( 16 ) with a finite PV. The inverse temperature 
P can be determined by fitting the distribution of the interac- 
tion energies that follow a Boltzmann law, i.e., the probability 
/?£ to find the system in a state with interaction energy 8 is 
/?£ oc e~^^. Further numerical analysis shows that the quantum 
state at long times not only shows this characteristic distribu- 
tion of the interaction energies, but also that other elementary 
relations known from thermal equilibrium states such as the 
fluctuation-dissipation theorem are satisfied 

It will be interesting to explore this further, not only the- 
oretically, but in particular experimentally. This will shed 
light on the question: Under which conditions can aspects of 
a closed many-body quantum system be understood by equi- 
librium thermodynamics? Experiments with Rydberg atoms 
in large lattices will be especially useful to probe and charac- 
terize the nature of the transition occurring in the vicinity of 

Vtrans • 
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4 Quantum state preparation 

The subject matter of this section is again the model ^ but 
with focus on the weakly interaction regime, i.e., V ^ Q.. 
We will see that in this regime the dynamics of the system 
takes place in constrained subspaces in which the quantum 
evolution is approximately governed by an exactly solvable 
Hamiltonian. The properties of the eigenenergies and the 
corresponding many-particle eigenstates will be analyzed and 
a scheme for their experimental realization will be outlined. 
Finding simple ways for creating entangled many-particle 
states is of importance for numerous applications: They serve 
as resource for precision quantum measurements^^, for mea- 
surement based quantum computing as well as for the cre- 
ation of single-photon light sources^D. We will conclude 
by discussing in particular the realization of a single-photon 
source based on a delocalized excitation within the Rydberg 
ring. 

4.1 Constrained dynamics 

In the weakly interacting regime the laser Hamiltonian ^ (in 
particular the term proportional to the Rabi frequency, since 
we assume Q,:^ |A|) dominates the dynamics. It is therefore 
convenient to choose a basis in which it is diagonal. This is 
achieved by the unitary transformation 



(17) 



and 



which brings G 
Hamiltonian ^ we obtain 



-Gj. When applied to our 



= U^HU = 



VL 



with 



xy 

Hi 
Hi 



L 

I 

^=1 L 



G_G\ 



= ft (■-■'?•) 



k=\ 
L 



(k) 



(18) 

(19) 
(20) 
,(21) 



where is the well-known x}^-model of a chain of spin 1/2 
particles in a transverse magnetic field. 

Let us now analyze the importance of the individual contri- 
butions to H\ As we can see in Fig. [9j the spectrum of the 
Hamiltonian decays into manifolds of states which are sep- 
arated by gaps whose width is approximately 2Q.. This is 

(k) 

caused by the dominant first term of //xy, i.e., Q'Y.k^z • The 



m+2 



m+1 




Fig. 9 Level structure in the regime Q.^V and |A| <C ^1. The 

spectrum splits into manifolds which can be labeled by the quantum 

(k) 

number m of the operator Y^k^z ■ For sufficiently large Q., the 
coupling between manifolds that is established only by Hi and H2 
can be neglected. The (constrained) dynamics inside the 
m-subspaces is then solely determined by //xy. 



(k) 

eigenstates of gI are - in terms of the (super)atom states 
given by 

1 



|±),= — //t[|G),±|/?),] 



with a^"^ = ± Thus, each of the manifolds that de- 
termine the coarse structure of the spectrum is spanned by a set 
of product states that have the same number of (super)atoms in 
the state | +) . In Fig. |9]we show the corresponding coarse level 
structure. There, we label each manifold by the eigenvalue of 
its states with respect to the operator YJkGl \ m, which basi- 
cally counts the difference between the number of sites that 
are in the |+) and the |— ) state. 

The second term of H^y conserves the total number of |+) 
(super)atoms. In other words, it couples only states that be- 
long to the same m-manifold and that are nearly degenerate. 
As a consequence, the strength of these intra-manifold cou- 
plings due to //xy is proportional to the interaction strength V. 
Conversely, Hi and H2 couple states that belong to manifolds 
with different number of (super)atoms in the state |+). In par- 
ticular. Hi and the last term of H2 flip one of the (super)atoms 
from 1+) to |— ) or viceversa. Thus, the coupled states belong 
to different manifolds with Am = ±1, energetically separated 
by The two first terms of H2 drive a similar process, flip- 
ping always two contiguous (super)atoms in the same state 

simultaneously, i.e., |++) \ ) or | ) |++). As a 

result, these terms connect states with eigenvalue m to those 
with m ± 2 and which possess a energetic separation of ap- 
proximately These selection rules are indicated in Fig. 
[9] The transition rates between m-manifolds corresponding to 
Hi and H2 can be estimated by second order perturbation the- 
ory to be of the order A^/Q and V^/Q., respectively. Hence, 
for sufficiently strong driving Cl^ V, their contribution can 
be neglected and the system's dynamics is constrained to the 
m-manifolds. As a consequence, the Hamiltonian that drives 



8| 



the intra-manifold dynamics, given by H^y, effectively drives 
the dynamics of the entire system in the considered parame- 
ter regime. This Hamiltonian is analytically solvable, and we 
thus have access to the actual spectrum and eigenstates of the 
system. 

4.2 Many-particle eigenstates 

We will now solve the Hamiltonian in order to determine 
the eigenstates and eigenenergies. This is done using the so- 
called Jordan- Wigner transformation and a subsequent Fourier 
transform^. 

The Jordan- Wigner transformation introduces the operators 



The ground state is given by 



|0) 



L 

n 

k=\ 



(25) 



and the excited states are formed by the successive application 
of the creation operators on the ground state. In particular, the 
singly and doubly excited states are of the form 



|7> 



l'7> 



«;|0). 



(26) 



(22) 



which are highly non-linear in the spin operators and obey the 
canonical fermionic algebra 



{c],Cj} = dij {cj,c]} = {Ci.Cj} 



0. 



After this transformation, the Hamiltonian (T9\ takes on the 
form 



xy 



2£1 



L 

E 

k=i 

^ ( t 

4 



clck- 



-c\cl 



respectively. Note that the notion 'singly and doubly excited' 
does not mean that one or two Rydberg atoms are excited. It 
means that one or two out of the L atoms on the ring are in the 
state 1+) while all others are in | — ). This means that the mean 
number of Rydberg atoms forming all the above-mentioned 
many particle states is L/2. 

4.3 Experimental preparation 

Our aim is now to find a scheme which allows to populate the 
excited many-body quantum eigenstates of the system in an 
experiment. In the same spirit as in Sec. |3j we assume that the 
initial condition is such that no Rydberg atom is present and 
the system is thus prepared in the state |init) =\{k\G)^. The 
proposed scheme consists of two steps, the preparation of the 



•1 



(23) 



ground state ( 25 ) of the Hamiltonian ( 24 ) and the subsequent 
excitation of the states ([26l). 



Hence, the Hamiltonian has been transformed into one which 
describes a chain of spinless fermions with nearest neighbor 
hopping. The last term of Hamiltonian ( [23] ) appears due to 
the periodic boundary conditions. It depends on the operator 



4.3.1 Preparation of the ground state of H^yi In the fol- 
lowing we will consider the simple situation in which there is 
only a single atom per lattice site. In this case the superatom 
state \R)^ becomes simply the single-atom state |r)^. We now 



which counts the total number of fermions, introduce, in addition to the single-atom ground state a 



which is also equivalent to the number of sites in the state | + ) . 
We thus have to explicitly distinguish between an odd or even 
number of fermions (which we refer to as parity). With this 
distinction being manifest, H^y reads 



(e/o) 
xy 



^=1 



second stable state |^)^. The latter will in practice be consti- 
tuted by a hyperfine state of the atomic ground state manifold 
different from | g ) ^ . 

Using this additional state we can prepare the ground state 



yL-l 

4 



k=\ ^ 



(25) starting from the initial state |init) by a sequence of two 
resonant laser pulses with Rabi frequencies 0^\ 2- We choose 
the first pulse to be resonant on the single atom transition 



V 



{c\ci^c\ci)j , 



\s) 



for a time Xi = 7l/ {2Q.i). Subsequently we irradi- 



^4 l^L^i 



for even (e) or odd (o) parity, respectively. 

We proceed by introducing the Fourier transformed opera- 



tors Tj];^ = 



Li 



exp (iank) c[ with the index p = o/q de- 



noting the parity and with = 2nn/L and = (2^ — l)7l/L. 
This finally leads to the diagonal representation of the Hamil- 
tonian: 



iffy = -La^Y. '^U'^P.n + cos < 



n=\ 



V 



(24) 



ate a strong laser (^^2 ^ ^) that resonantly couples |^)^ |r)^ 
for a time X2 = ^/^2- This amounts to the sequence 

linit) = ni^)^ (27) 

k 

^ X{^[\s)k + i\s)k\ 

^ X{^[\s)k-\r)k\=Y[\-)k = \^) 
k v2 

and thus results in the desired preparation of |0). 
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4.3.2 Excitation of many-body quantum states from 

|0): Let us show now how to address the single-fermion and 
two-fermion states from |0). To this end we make use of the 
Hamiltonian H\ ([20|) which emerged from the laser detuning 



after the appHcation of the unitary transformation ([T7]). As we 
discussed earher, this term drives transitions between mani- 
folds with Am = ±1, (see Fig. |9]) and such a transition would 
exactly lead to the excitation of the desired many-particle 
states. However, we have also estimated earlier that these tran- 
sitions have small probability since they are suppressed by a 
factor - with |A| <C O^. 

To overcome this problem we introduce an oscillating de- 
tuning of the form A(^) = Aqsc cos (cOaO- tunes cOa 
such that it coincides with the gap between two given states, 
this detuning acts effectively as a laser that couples them res- 
onantly (see Fig. [10^). The matrix element corresponding to 
this transition (within the rotating-wave approximation) and 
thus the transition rate is given by 



(final|ifi|0) = ^(final|X;cji')|0). 

^ k=\ 



(28) 



It turns out that, due to the symmetry of the Hamiltonian and 
the operator H\, this matrix element is non-zero only for a 
single final state given by 



whose energy is 



V 

r 



(29) 



(30) 



with = —L (^1 — ^ ) being the ground state energy. This 
state is a spin wave or, in other words, a superatom that ex- 
tends over the entire lattice. 

Once this spin wave is excited, one can use the oscillating 
detuning to reach states that carry two excitations. Again, the 
selection rules restrict the number of accessible states and only 
the transitions to 



(31) 



iL 



sm 



k>k' 



2% 



afaf)|0) 



are permitted. The energy of these states is given by 

"271, 



= £o+4n+ycos 



-{p-\/2) 



(32) 



A schematics of the envisioned excitation path is depicted in 



I2p> b 
|1> 

|0> 



CO, 



12^ 
|1> 

|0> 



Fig. 10 Sketch of the excitation of the single-fermion and 
two-fermion states by means of an oscillating (radiofrequency) 
detuning using a not too large value of 12. a: In a first step, the 
population is transferred by a 7i-pulse to the single-fermion state by 
tuning the frequency of the detuning on resonance with the gap 
COa = coi = ^1 — ^0- b» A second 7C-pulse with coa tuned to match 
CO2 = E2p — El addresses the corresponding |2p) state, bearing in 
mind that V ^ Aqsc in this step. 



scheme. Firstly, in order to avoid accidental resonances during 
the excitation process, one can exploit the second-order level 
shifts that are caused by H2 in the regime in which the ratio 



is not too small. This is sketched in Fig. 10 i, where the 
gap between 1 1) and any of the \2p) levels, i.e., CO2 = — E\ , 
becomes increasingly different from the energy gap between 
|0) and i.e., COi, for sufficiently large V /Q.. Secondly, the 
separation between neighboring \2p) states is of the order of 
y and the Rabi frequency of the transition is proportional to 
Aqsc- As a consequence, to populate only a single level of 
the doubly-excited manifold, the parameters have to accom- 
plish that y ^ Aqsc and, at the same time, Aqsc has to be large 
enough in order to perform the transfer in a time interval that 
is much shorter than the lifetime of the Rydberg state. 

4.4 A single photon source 

A collective atomic excitation stored in the Rydberg ring can 
be converted into photonsS^^^Sl xhe particular features of the 
atomic state are reflected in the emission properties, i.e., the 
angular distribution of the emitted radiation. We will illustrate 
this using the single excitation state ([29]) as an example. 

We start by mapping the excited state which is encoded 
in the superposition states |+)^ and |— )^ to stable states of 
the atomic hyperfine groundstate manifold. This is necessary 
since such states are much longer lived than the Rydberg states 
which decay due to black-body radiation or spontaneous emis- 
sion. Reverting the sequence ([27]) - which led to the prepara- 
tion of the ground state |0) of H^y - performs the mapping 



-)k 



\8)k 
i\s)k- 



(33) 
(34) 



Fig. 10 There are some limitations to bear in mind in this 



In order to convert the collective atomic excitation stored 
in the Rydberg ring into photons we make use of the level 
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A — li>k 

C\_e'''J \ |k|,e,(t) 
\ 

\ 

\ 

— |gX 



Fig. 11 a: Level scheme for the photon generation scheme. The 
delocaUzed atomic excitation is stored in the local stable states 
and 1^)^ which is achieved by a suitable mapping sequence 
converting states |±)^ into ground states. The state is coupled 
by a laser with wavevector ki and Rabi frequency Q.l off-resonantly 
(detuning A^) to the intermediate state |/)^. The photon is then 
emitted on the transition \g)]^. b: Outcoupling scheme. We 

consider a situation in which the atomic transition dipole moments 
are parallel to the z-axis (small arrows) and the wavevector of the 
laser is k/^ = I^lIg^- Photons are then emitted into modes with a 
certain angular distribution, parameterized by the angles 9 and (|). 



scheme depicted in Fig. 11 i^U. Here a laser with wavevec- 
tor couples the state \g)j^ off-resonantly to an intermediate 
state |/)^. The Rabi frequency and the detuning of this tran- 
sition are given by Q.^ and A^, respectively. Photons are then 
emitted on the transition \g)]^. In our scheme we do not 

consider decay from |/)^ back to l^-)^, whic h can be ensured by 
an appropriate choice of atomic 

levelsSHSl. 

We consider now times much longer than the lifetime of 
the intermediate state x = F"^, where F is the corresponding 
decay rate. One can show that in the above described scheme 
a single atomic excitation is mapped to a single photon state 
according to 



^L^J-qv^qvlvac) 



(35) 



where a^^ creates a photon with momentum q and polariza- 
tion V and I vac) is the photon vacuum. The coefficients gjqy 
contain information about the coupling strength of the excita- 
tion on the 7-th lattice site to the photon mode characterized 
by q and v, and are given by 



Sjqv{t) : 



-iqvy f^J 



/CO - Dk 



(36) 



Here, CO = |q|/c, ry = R (cos(|)y, sin4)Y,0) with (^y = ^(y- 1) 
denotes the position of the atoms and the coefficient reads 



qv 



Al 



CO 



dgi e, 



qV7 



(37) 



where V is the quantization volume, 8o the vacuum permitiv- 
ity, and dgi is the dipole operator of the \g)j^^ transition. 



!A4 and Dk are the eigenvectors and eigenvalues of the oper- 
ator that governs the atomic dynamics and which we call /. 
It depends on the orientation of the atomic transition dipole 
moments and the quantity k^R, where aL/2Kis the radius 
of the ring. We consider the particularly simple situation in 
which the transition dipole moments are aligned with the z- 
axis (see Fig. [TTJ)). In this case, the matrix representation of / 
is a circulant complex symmetric matrix and its eigenvalues 
and eigenfunctions are given by 



where J\n 



w-^ -In 



(38) 



and 



3F 


"cOSKi^ 


sinKi^ 


1 sinKi^" 


Y 


[ < 






3F 


"sinKi„ 


COSKi^ 


COSKi^" 


Y 









withKi„=^L/?|ri-f^P. 

To study the angular distribution of the emitted photon we 
study the intensity of photons per solid angle that is defined 
through 



/(e,(^) 



V 



{2ncy- 



[ y (^qv)co^^o). (39) 
^0 V 



In the particular case of the spin wave given in (29) and the 
direction of being parallel to the z-axis, one can prove that 
the expression of the angular intensity yields 



/(e,(^) 



3F sin^e 



47lLZ)i 



y=l 



-ikiRq-fy 



(40) 



with q = (sin cos (|), sin sin (|), cos 0). For sufficiently small 
values of kia, i.e., kia ^ 7r/2, the sum in the previous expres- 
sion can be substituted by a Bessel function and the intensity 
becomes 



7(0,^); 



3FL sin^0 
~47Z Di^D 



-J^ihRsind). 



(41) 



In Fig. 12 we show the intensity for kia = (0.25, 1.36, 6.0) 
and L = 40. The figure shows strikingly how the distance 
between the atoms affects the profile of the radiation. For 
kia = 0.25 the atoms are so close together that the atomic exci- 
tation (spin wave) acts as a single degree of freedom that cou- 
ples to the radiation field This results in an almost spher- 
ical intensity profile which is modulated by the dipole radia- 
tion pattern. In the intermediate case kia = 1.36 the coupling 
to the radiation field is still dominated by collective effects. 
Here, the photon emission is strongly peaked along a polar 
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kLa=0.25 



Ka=1.36 



kLa=6.0 




Fig. 12 Angular intensity distribution for a ring of L sites and 
various values of kia. The density plots in the lower row show the 
same data as the three-dimensional plots. The red arrow indicates 
the direction of \ii which is chosen to be parallel to the z-axis. For 
kid < 7i/2 the intensity profile is created through the collective 
coupling of the atoms to the radiation field. Here no dependence on 
the azimuthal angle is visible. For larger kia the atoms couple 
individually to the electromagnetic field which gives rise to a large 
number of peaks. 



angle ^ 7r/3. The position of this peak can be controlled 
by tuning the parameter kia. In the previous two cases, the 
intensity profiles show no dependence on the azimuthal an- 



gle, as we expected from expression (41 ). In the third case, 
kia = 6.0, the equation (41 ) is no longer valid, and the atoms 
can be approximately regarded as independent so that they are 
coupled individually to the radiation field. This produces a 
large number L of peaks, as a result of the interference of the 
L sites. 

In addition to these examples the Rydberg ring also allows 
to generate single photons with a strong directionality, i.e. the 
photon emission takes place only into a very small solid angle. 
Moreover, by using many-body states that contain more than 
a single excitation, the creation of correlated photon pairs or 
triples can be achieved 

5 Conclusions and Outlook 

In this paper we could merely offer a slight glimpse on the 
potential the highly excited atoms have to offer for the study 
and the understanding of quantum many-body physics. We 
saw that already the very simple Rydberg ring allows to tackle 
such fundamental questions as the thermalization of a closed 
quantum systems. Moreover, we demonstrated that the Ry- 
dberg ring can be used to create single photon sources with 
particular emission characteristics. 



More interesting features are expected to emerge in higher 
dimensional setups, particularly in two dimensions where 
many-body systems possess a particularly rich phase struc- 
ture. While certainly some of the concepts that were presented 
here are also applicable in higher dimensions there will also be 
fundamental differences. E.g., in the limit of strong laser driv- 
ing the system will no longer be analytically solvable as the 
Jordan- Wigner transformation leads only to a solvable model 
in one dimension. 

The current experimental situation is such that Rydberg 
atoms are excited from an atomic gas in which the atoms are 
randomly distributed, and the experimental realization of the 
Rydberg ring is certainly several steps further down the road. 
Significant efforts are, however, undertaken to study the Ryd- 
berg excitation dynamics in structured environments and first 
exciting results have shown the feasibility of these undertak- 



ings 
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